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Recently measured lusion cross-sections for the neutron-rich system 20 O+ 12 C are compared to 
dynamic, microscopic calculations using time-dependent density functional theory. The calculations 
are carried out on a three-dimensional lattice and performed both with and without a constraint on 
the density. The method has no adjustable parameters, and its only input is the Skyrme effective 
NN interaction. While the microscopic DC-TDHF calculations lie closer to the experimental data 
than standard fusion systematics they underpredict the experimental data significantly. 
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The outer crust of an accreting neutron star provides 
a unique environment in which nuclear reactions can oc- 
cur. It has been proposed that the fusion of two neutron- 
rich light nuclei in the outer crust could provide a heat 
source to ignite thermonuclear fusion of 12 C + 12 C and 
produce a signature X-ray superburst 1. To date, how- 
ever, a limited amount is known either experimentally or 
theoretically about the fusion of neutron-rich nuclei. Pi- 
oneering experiments with heavy nuclei indicate that the 
fusion below the barrier may be enhanced [2j [3] . Such an 
enhancement has recently been associated with the im- 
portance of neutron transfer channels which effectively 
lowers the fusion barrier [4 . In the case of fusion of two 
neutron-rich light nuclei (Z<20), even less is known. In 
principle, this is the most promising domain as neutron- 
rich nuclei up to the drip line can be experimentally pro- 
duced. 

Recent experimental measurement of near-barrier fu- 
sion in the system 20 O + 12 C [5 suggests that the fusion 
cross-section is enhanced relative to the predictions of 
the Bass model [6] . As the empirical Bass model is based 
upon the systematics of known fusion cross-sections near 
/3-stability, it does not include the increased importance 
of neutron transfer channels for neutron-rich nuclei. The 
aim of this paper is to directly compare the experimen- 
tal results with a microscopic approach, namely the time 
dependent Hartree-Fock (TDHF) theory. 

The experiment was performed at the SPIRAL1 facil- 
ity at the GANIL accelerator complex in Caen, France. 
An 20 O beam with an intensity of 1 — 2 x 10 4 p/s im- 
pinged on a 100 jug/cm 2 thick 12 C target. The energy 
of the beam on target was varied between 1 MeV/A and 
2 MeV/A in order to measure the fusion excitation func- 
tion. Experimental details have been previously pub- 
lished [5 j and are summarized here only for completeness. 
Nuclei produced by fusion subsequently de-excite via 
evaporation of neutrons and light charged particles (Z<2) 
forming evaporation residues. These residues were de- 
tected in two segmented silicon detectors located down- 
stream of the target and identified by measuring both 



their energy and time-of- flight. The annular detectors 
spanned the angular range 3.54° < 9i a b < 21.8°. Due to 
the presence of a large atomic background in the experi- 
ment, a coincidence between an emitted charged particle 
and the evaporation residue was necessary to distinguish 
fusion reactions. Statistical model calculations with a 
Hauser-Feshbach model, evapOR, indicate that, depend- 
ing on the excitation energy of the compound nucleus 
formed, approximately 15-25 % of fusion reactions de- 
excite via emission of at least one charged particle. 

The time-dependent Hartree-Fock (TDHF) theory pro- 
vides a useful foundation for a fully microscopic many- 
body theory of large amplitude collective motion. It is 
therefore well suited to describing deep-inelastic and fu- 
sion reactions [7, 8 . Only in recent years has it become 
feasible to perform TDHF calculations on a 3D Cartesian 
grid without any symmetry restrictions and with much 
more accurate numerical methods [8 13 . In addition, 
the quality of effective interactions has been substantially 
improved [T4T - H7] . TDHF theory predicts an energy den- 
sity functional which is determined by the given effective 
NN interaction. One may therefore view TDHF as a 
special case of a time-dependent density functional the- 
ory (TDDFT), a concept used in many areas of nuclear 
physics, condensed-matter physics, and chemistry. 

Over the past several years, the Density Constrained 
Time-Dependent Hartree-Fock (DC-TDHF) method for 
calculating heavy-ion potentials [TH] was utilized to cal- 
culate fusion cross-sections. We have applied this method 
to calculate fusion and capture cross-sections above and 
below the barrier to about 20 systems to date, exam- 
ples of which can be found in Refs. [19 23 . Recently, we 
have also investigated sub-barrier fusion between nuclei 
that occur in the neutron star crust [4]. In all cases, we 
have found good agreement between the measured fusion 
cross-sections and the DC-TDHF results. This agree- 
ment is rather remarkable given the fact that the only 
input in DC-TDHF is the Skyrme effective N-N interac- 
tion, and there are no adjustable parameters. 

The TDHF equations for the single-particle wave func- 



tions 



d 



H{<f>^})Mr,t) = ih- ( f )x (r,t) (X=1,...,A), (1) 

can be derived from a variational principle [7]. In the 
present TDHF calculations we use the Skyrme SLy4 in- 
teraction [2] for the nucleons including all of the time- 
odd terms in the mean- field Hamiltonian |10| . The nu- 
merical calculations are carried out on a 3D Cartesian lat- 
tice. For the calculations shown in this work, the lattice 
spans 40 fm along the collision axis and 24 — 30 fm in the 
other two directions, depending on the impact parameter. 
We first generate very accurate static HF wave functions 
for the two nuclei on the 3D grid. In the second step, we 
apply a boost operator to the single-particle wave func- 
tions. The time-propagation is carried out using a Tay- 
lor series expansion (up to orders 10 — 12) of the unitary 
mean-field propagator, with a time step At = 0.4 fm/c. 

Presented in Fig. [T] is a contour plot of the mass den- 
sity during a collision which clearly shows the formation 
of a neck between the two fragments. This density dis- 
tribution, shown here for 20 O+ 12 C at £ c . m . = 9.5 MeV, 
is representative of collisions for similar systems. As the 
collision proceeds in the TDHF calculation, transport of 
protons and neutrons between the two nuclei can be fol- 
lowed within the theory. For larger impact parameters 
the larger angular momentum of the system leads the 
two nuclei to separate and a deep-inelastic reaction oc- 
curs. For smaller impact parameters the disrupting in- 
fluence of angular momentum and Coulomb repulsion is 
insufficient to overcome the nuclear attraction and fusion 
results. By examining the density distribution as the two 
nuclei fuse into one within the calculation, one clearly ob- 
serves the occurrence of a damped dipole resonance and 
surface waves. Deep inelastic and fusion reactions are the 
dominant reaction channels in this energy domain. Dis- 
tinguishing between these two types of reactions is real- 
ized by examining the density distribution as a function 
of time and observing whether one or two large fragments 
result from the collision. 

In the absence of a true quantum many-body theory of 
barrier tunneling, all current sub-barrier fusion calcula- 
tions assume the existence of an ion-ion potential V(R) 
which depends on the internuclear distance R. Most of 
the theoretical fusion studies are carried out with the 
coupled-channels (CC) method [241127] in which one uses 
empirical ion-ion potentials (typically Woods-Saxon po- 
tentials, or double-folding potentials with frozen nuclear 
densities). In contrast to the use of these empirical poten- 
tials we have adopted a microscopic approach to extract 
heavy- ion interaction potentials V{R) from the TDHF 
time-evolution of the dinuclear system which describes 
the dynamics of the underlying nuclear shell structure. 

In the DC-TDHF approach [HJ, the TDHF time- 
evolution proceeds uninhibitcdly. At certain times t or, 
equivalently, at certain internuclear distances R(t) the 
instantaneous TDHF density is used to perform a static 
Hartree-Fock energy minimization while constraining the 
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FIG. 1. (Color online) Unrestricted TDHF+BCS calculation 
for 20 O+ 12 C at £" c . m . = 9.5 MeV and impact parameter 6 = 
2.5 fm. Shown are mass density contours shortly after a neck 
has formed between the two fragments. 



proton and neutron densities to be equal to the instan- 
taneous TDHF densities. This means that we allow the 
single-particle wave functions to rearrange themselves in 
such a way that the total energy is minimized, subject 
to the TDHF density constraint. In a typical DC-TDHF 
run, we utilize a few thousand time steps, and the den- 
sity constraint is applied every 10 — 20 time steps. We 
refer to the minimized energy as the "density constrained 
energy" Evc{R)- The ion-ion interaction potential V(R) 
is obtained by subtracting the constant binding energies 
E Al and E^ 2 of the two individual nuclei 



V(R) = E DC (R)-E Al -E A2 



(2) 



In direct TDHF calculations the fusion cross-section 
is calculated by determining the maximum impact pa- 
rameter for which fusion occurs and applying the sharp 
cut-off approximation. For example, in the case of the 
reaction 20 O+ 12 C at E c . m . = 9.5 MeV wc find that im- 
pact parameters b < b max — 4.075 fm result in fusion, 
while impact parameters b > b max lead to deep-inelastic 
reactions. Using the sharp cut-off model, the fusion cross- 
section is given by CTf us = ^b 2 max = 52.2 fm 2 = 522 mb. In 
contrast, in DC-TDHF method fusion cross-sections are 
obtained by integrating the Schrodinger equation for the 
potential V(R) with a coordinate-dependent mass [2Tj . 

We begin the comparison of the microscopic calcula- 
tions with experimental data by examining the well stud- 
ied system 16 0+ 12 C. Shown in Fig. [2] are four sets of ex- 
perimental data for the total fusion cross-section along 
with the corresponding microscopic DC-TDHF calcula- 
tions. The experimental techniques used to determine 
the fusion cross-section range from gamma spectroscopy 
of the charged particle channels (Christensen and Cu- 
jec) to direct measurement of the evaporation residues 
(Eyal and this work) . It should be stressed that the data 
represented by the red filled circles (this work) utilized 
the same experimental technique and setup as for the 
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FIG. 2. (Color online) Comparison of the experimentally mea- 
sured fusion excitation function with the predictions of the 
DC-TDHF method for 16 0+ 12 C. Data from [213H5U1 shown. 



20 O+ 12 C results subsequently presented. For energies 
E c .m. < 9.5 MeV all the experimental data are in agree- 
ment. At higher energies however, the cross-sections 
measured by Christensen et al. slightly exceed that of 
the three other datasets. We have no explanation at 
present for the larger cross-sections measured by Chris- 
tensen et al. Since the cross-section measurements of 
Eyal, Cujec, and the present work are all in good agree- 
ment we take these cross-sections to accurately represent 
the true fusion cross-section. It is interesting to note that 
the DC-TDHF calculations also slightly exceed the mea- 
sured cross-sections for E cm > 7 MeV. At the largest 
energies measured this excess is of the order of 20 %. In- 
deed, the agreement of DC-TDHF results with the data 
of Christensen et al. at the highest energies was some- 
what surprising since TDHF dynamics for light heavy- 
ions at these energies do not properly account for various 
breakup channels present for these systems, and results in 
a fusion-like composite system with long-time collective 
oscillations. In coupled-channel calculations this discrep- 
ancy is cured by introducing a small imaginary potential 
in the vicinity of the potential minimum |27j . Having 
established the degree of confidence through the com- 
parison of the fusion cross-sections in 16 0+ 12 C, we then 
calculated fusion in 20 O+ 12 C. 

While for isolated 16 and 12 C nuclei the Hartree-Fock 
(HF) ground state is found to be spherical in agreement 
with both theory and experiment, for 20 O the HF cal- 
culations predict a prolate quadrupole deformation of 
02 = 0.25. This deformation is in disagreement with self- 
consistent mean field calculations with pairing (Skyrme- 
HFB) [31j which predict a spherical nucleus. Moreover, 



the measured energy level spectrum [32] also shows this 
nucleus to be spherical. In addition, a measurement of 
the magnetic moment of 20 O 03] also indicates its spher- 
ical nature. We attribute this prolate deformation pre- 
dicted by the HF calculations to the lack of pairing in 
the method. One unfortunate consequence of the DC- 
TDHF approach in calculating the fusion cross-section 
lies in the treatment of pairing during the collision pro- 
cess. In unrestricted TDHF calculations, the BCS occu- 
pation numbers can be kept frozen during the collision 
to have correct initial states. This approximation cannot 
be utilized in the DC-TDHF method because the static 
HF solution coupled with a constraint on the instanta- 
neous TDHF density for the combined system requires 
the reevaluation of the occupation numbers for the low- 
est energy solution. Consequently, calculations with the 
DC-TDHF method do not include pairing. On qualita- 
tive grounds, it can be argued that this omission should 
result in a slightly larger prediction of the fusion cross- 
section. The reason is that pairing results in a spherical 
20 O nucleus, and the fusion barrier for a spherical nu- 
cleus is higher than the lowest barrier for a deformed 
nucleus. In order to calculate the fusion cross-section 
for this system within DC-TDHF method we therefore 
take an average of all initial orientation angles /3 of the 
deformed 20 O nucleus, where f3 is defined as the angle be- 
tween the internuclear distance vector and the symmetry 
axis of the deformed nucleus. 

As the collision occurs, using TDHF dynamics, it is 
possible to compute the corresponding coordinate depen- 
dent mass parameter M(R) [21]. At large distance R, the 
mass M(R) is equal to the reduced mass fj, of the sys- 
tem. At smaller distances, when the nuclei overlap, the 
mass parameter generally increases. In order to calculate 
the fusion cross-section more easily, one can replace the 
coordinate-dependent mass M{R) and the original poten- 
tial V(R) with the constant mass fj, and the "transformed 
potential" U(R), using a scale transformation [21]. In 
Fig. k| we display the transformed potentials U(R) for 
initial orientation angles /3 = 0°, 10°, ...,90° of 20 O. For 
sufficiently large separation between the two nuclei, R 
>9, fm all the transformed heavy-ion potentials are the 
same regardless of the orientation of the two nuclei. With 
decreasing distance as the two nuclei come into contact, 
the heavy-ion potentials differ. For a nucleus with pro- 
late deformation, the orientation angle /3 = 0° leads to 
the lowest potential barrier. This reduction in the poten- 
tial occurs because the distance between the nuclear sur- 
faces is minimized for this orientation which is important 
because of the short-range nature of the strong interac- 
tion. The Coulomb interaction is accurately calculated 
by solving the 3D Poisson equation numerically during 
the collision. We observe that with increasing orienta- 
tion angle, the barrier height increases and the barrier 
position is shifted to smaller distances. 

Displayed in the inset of Fig. [3] is the fusion excitation 
function associated with the orientations of /3=0° and 90° 
as well as the angle averaged result. As one might quali- 
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FIG. 3. (Color online) Transformed heavy-ion potentials 
U(R), corresponding to the constant reduced mass /j,. The 
potentials have been calculated for a series of initial orienta- 
tion angles j3 of the deformed 20 O nucleus (no pairing). Shown 
in the inset is the total predicted cross-section as a function of 
E c .ta. for /3=0°, /3—90°, as well as the cross-section averaged 
over p. 



tatively expect, based upon the heavy-ion potentials, for 
a given E c . m ., the orientation /3 = 0° with the lowest 
potential barrier corresponds to the largest cross-section. 

In Fig. [4] we compare the microscopic calculations us- 
ing the DC-TDHF method to the experimental data. The 
experimental data at £ c . m = 7.35, 9.29, and 15.24 MeV 
are shown as the filled circles. Due to the atomic back- 
ground previously mentioned only the fraction of the 
fusion cross-section associated with subsequent charged 
particle emission was experimentally measured. In or- 
der to compare the microscopic calculations with the ex- 
perimental data, we therefore calculated the fraction of 
compound nuclei produced at each incident energy that 
de-excite via emission of at least one charged particle. To 
calculate the de-excitation of the fused nuclei we utilized 
a Hauser-Feshbach statistical model, evapOR. The the- 
oretically predicted fusion cross-section associated with 
subsequent charged particle channels is depicted by the 
dashed line. This predicted cross-section clearly un- 
derpredicts the experimentally measured cross-section. 
At the highest energy, E cm =15.24 MeV, the predicted 
cross-section is 60 % of the experimentally measured one 
while at the lowest energy the predicted cross-section is 
substantially lower, only 30 % of the experimental value. 
For reference, we also present, as a dot-dash line, the fu- 
sion cross-section predicted by the Bass systematics that 
is associated with charged particle emission. This cross- 
section is less than that of the DC-TDHF method most 
likely reflecting the influence of neutron transfer chan- 
nels in aiding the fusion process. It should be appreci- 



ated that not only does the DC-TDHF method predict 
a larger cross-section at all energies as compared to the 
Bass systematics, but this increase grows with decreas- 
ing -E c .m.- This result suggests that the neutron transfer 
becomes more important in the sub-barrier domain. The 
result that the DC-TDHF method under/predicts the ex- 
perimental data by a significant amount is noteworthy. 
Moreover, it should be noted that the lack of pairing 
in the DC-TDHF method and the resulting deformation 
of the 20 O, as previously discussed, acts to increase the 
predicted cross-section implying that the discrepancy be- 
tween theoretical prediction and experimental data is at 
least as large as that evident in Fig. |4] 
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FIG. 4. (Color online) Comparison of the angle-averaged 
DC-TDHF fusion cross-sections with experimental data for 

20 O+ 12 C. 

In order to assess the impact of pairing on the mea- 
sured cross-sections more quantitatively we have per- 
formed unrestricted TDHF calculations (no barrier tun- 
neling) at energies above the barrier. These calculations 
were initialized with BCS/Lipkin-Nogami pairing for 20 O 
which resulted in a spherical nucleus, consistent with the 
experimental data. During the collision of the 20 O with 
the 12 C the BCS occupation numbers are kept frozen. 
The results of these calculations are presented as the 
open triangles in Fig. [4] A slight reduction in the total 
cross-section is evident. This reduction is of the order of 
5-20 % with the largest reduction for the lowest energy 
point calculated. From these calculations one can infer 
that the inclusion of pairing in the DC-TDHF calcula- 
tions should result in a slight reduction of the predicted 
cross-section, thus increasing the discrepancy with the 
experimental data as anticipated. 

It is tantalizing to speculate about possible reasons for 
the discrepancy between the theoretical predictions and 



the experimental data. In general terms, the discrepancy 
between the experimental data for the charged particle 
channels in Fig. 4] (solid points) and the corresponding 
DC-TDHF calculations can be thought of as originating 
either from an underprediction of the total fusion cross- 
section or from an underestimation of the relative impor- 
tance of the charged particle decay in the de-excitation 
of the compound nucleus. Either or both of these sources 
could explain the underprediction of the cross-section. It 
is therefore important to not only measure the total fu- 
sion cross-section but also the cross-section for individual 
decay channels. Furthermore, since the underprediction 
exists for energies well above the barrier, E ss 9.5 MeV, 
it is not simply a question of enhanced tunneling. More- 
over, beyond the overall underprediction of the micro- 
scopic method, one observes that the experimental data 



manifests a slower fall-off for the measured cross-section 
with decreasing incident energy as compared to the DC- 
TDHF calculations. This experimentally determined en- 
ergy dependence is highly provocative and it remains to 
be seen whether the total fusion cross-section also ex- 
hibits this slower fall-off. As this fall-off is intimately 
related to the neutron transfer channels, near and sub- 
barrier fusion of neutron- rich nuclei provides direct access 
to the extent of the neutron density distribution. 
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